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Chapter 1 



Fundamental Concepts 



1.1 Introduction 

Population transfer in atoms and molecules is one of the most intriguing problems of quantum optics. 
To carry out successful multilevel excitations is an important task for theoretical and practical 
purposes. Thus, figuring out different ways of improving the efficiency of those methods used to 
transfer population between multilevel systems has become an important subject of research during 
the past years. 

In this work, we examine numerically an atom + cavity system with two of the most important 
methods for transferring population between atomic levels: the 7r-pulses method and the adiabatic 
passage scheme. In particular, we explore the possibility of using the nonadiabatic characteristics of 
the second method to achieve very high transfer efficiencies without using large values of the pulse 
amplitudes or interaction times. 

The TT-pulse method takes advantage of the Rabi population oscillations that characterize coher- 
ent excitation. By adjusting the laser intensity and the pulse duration so that the time integral of 
the Rabi frequency (the pulse area) has the value tt, it is possible to produce complete population 
transfer between two states [SBK+92]. If the system undergoes no spontaneous emission, transfer 
efficiencies of 100% are possible to obtain with this method. However, a very accurate control of the 
pulse area is required, which is a really disadvantage. 

The adiabatic passage scheme is an alternative for transferring population between selected quan- 
tum states. The population can be transferred between two states by driving the system sufficiently 
slowly with the appropriate external fields, with the population remaining in an approximate energy 
eigenstate. This method is quite insensitive to changes in parameters like the laser pulse shapes, 
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intensity, and frequency modulation, as long as certain easily controllable experimental conditions 
are satisfied. It is important to note that the adiabatic passage method tends to minimize the pop- 
ulation of the intermediate level 2 in a three level system. Normally this level undergoes radiative 
decay. Therefore, the effects of spontaneous emission are largely reduced. For the adiabatic following 
solution to be valid, the detuning from the two photon resonance should be small compared to the 
Rabi frequency ilo- 

The dissertation opens with a mathematical and physical introduction of two- and three-level 
systems. Concepts like Bloch equations, rotating-wave approximation, the 7r-pulse method, and the 
adiabatic following are introduced here. This chapter is based mainly in the Allen and Eberly [AE75] , 
and Scully and Zubairy [SZ97] books. 

In the second chapter, we explore numerically an atom + cavity system with four atomic levels 
and a quantized coupling field. Adiabatic passage and 7r-pulse methods were used to drive the system 
classically. We find that for a very particular set of parameters, like Rabi frequency, pulse width, 
and time delay, we get very large transfer probabilities. 

In chapter three, we introduce a simple analytical model based in nonlinear differential equations, 
that can help to understand why for some nonadiabatic processes we still achieve very high transfer 
efficiencies. We found that the nonlinear system can be converted into a system of equations relating 
the angular coordinates of the state vector in the Hilbert space. In this way, we can explore how the 
state vector follows the adiabatic states of the Hamiltonian governing the evolution of the system. 

In the last chapter we consider the problem of coherence transfer of ground state levels, between 
two atoms inside a quantum microcavity We proved that adiabatic methods can be used for trans- 
ferring such coherence. In addition, we examined the system by using an alternative model based 
on symmetric and antisymmetric eigenstates, finding a connection between the original system and 
electromagnetic induced transparency and 27r-processes. 

1.2 The two-level atom 

Consider the solution to the Schrodinger equation for a two-level atom interacting with a classical 
coherent driving field. The state of the system is described in terms of the vector |^(i)), which 
obeys the time-dependent Schrodinger equation 

zfi||*(t))=i7|*(i)), (I.l) 



with the Hamiltonian operator H given by 



H = Ho-f^-'E{ro,t). 



(1.2) 



Here -ffo is the unperturbed Hamiltonian, fi is the atom's dipole moment operator, and E(ro,i) is 
the electric field operator evaluated at the position of the dipole. 

We assume that the applied electric field is quasi-monochromatic with a frequency nearly coinci- 
dent with the transition frequency connecting the atomic ground state \a) and some other level \b), 
as shown in Fig. 1.1. 
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Figure 1.1: The two-level atom coupled by a near- resonant transition. 

Because of the effect of the interaction is to mix states \a) and \b), we are only concerned with electric 
dipole transitions between these two levels. Thus the state vector of the system in the presence of 
the applied field can be written as 



\^{t)) = C„(t)e-'""*|a) + Cb(i)e-*"^*|6), 



(1.3) 



where Ca{t) and Ch(t) represent the probability amplitudes that at time t the atom is in state \a) 
or state |6), respectively. Assuming that these states are eigenvectors of the Hamiltonian Hq with 
eigenvalues huja and hujb, the matrix elements of the atomic operators can be written as 



{a\Ho\a) = hwa 
{b\Ho\a) - 



{a\H„\b)^0 
{b\Ho\b) ^ hiOb, 



(1.4) 



and 



{a\n\a) = {a-lfJ-lb) = fiab 

(6|M|a) = K, (6|m|6) = 0. 



(1.5) 



We note that there are no ofF-diagonal elements of Ho, because \a) and \b) are considered to be 
orthonormal eigenvectors of _ffo; and there are no diagonal elements of /x, because it is a vector 
operator and then has odd parity (we assume that \a) and |6) have definite parity). 
In general, the dipole matrix elements are complex vectors that can be expressed as 






(1.6) 



where iir and jjbi are real vectors. Therefore, we can represent the Hermitian operator fi by the 
two-dimensional ofF-diagonal matrix: 



M 



/J-r + ^/J-i 

fir — ifJ'i 



By introducing the two-dimensional Pauli matrix operators 

V 



0-1 



1 

1 



(^2 = 



-i 
i 



0-3 = 



1 



-1 



the unperturbed Hamiltonian and the atom's dipole moment operators can be written as 



-^0 = ^{^b+ ^a)I + ^{^b- ^a)<73, 



/X = /XrO'l — (J"i<y2- 



(1.7) 



(1.8) 



(1.9) 
(1.10) 



In this way, the Hamiltonian of the system takes the form 



H = -{uJb +Wa)/+ -{oJb -^a)a-3 - ifJ'r ' E)cri + (/Xj •E)cr2. 



(1.11) 



Here / represents the 2x2 identity operator. 

In the Hcisenberg picture, the equation of motion for an operator that does not depend explicitly 



on time is given by 



ihj^A^[A,H]. (1.12) 



Then, by plugging the Pauli matrix operators and Eq. (1.11) into the Heisenberg equation (1.12), 
we may obtain 

di{t) = -woff2(i) + J^[^^^■ E(t)]CT3(0, (l-13a) 

d2{t) = LJo<yi{t) + j:[f^r- E{t)]<J3{t), (1.13b) 

'^sW = -7:[Atr • E(t)](J2(t) -T[^^^■ E(t)]ai(i), (1.13c) 



where luq = (uh — LJa)/fi represents the atomic transition frequency, and E(i) has been taken in the 
Heisenberg picture. 

As we see from Eqs. (1.13), the operator nature of the atom and field variables makes the 
system very difficult to solve and no general solutions are known. Moreover, if the operator Maxwell 
equations that govern the electric field are included, the system becomes even more complicated. 
One way to overcome this problem is by considering an alternative system of equations for which 
the quantum correlations between field and atom can safely be ignored [AE75]. As a result, the 
expectation value of any product of two operators of the form E{t)ai{t) can be expressed as the 
product of the individual expectation values: 

(E(i)a.(t)) = (E(i))(a,(t)). (1.14) 

In this way, we can reformulate the semi-classical radiation theory of two-level atoms by applying 
this factorization systematically to Eqs. (1.13). Now, according to the following notation 

s^{t) = {<j,{t)), z = l,2,3 (1.15) 

E(i) EE (E(t)), (1.16) 

the set of equations (1.13), that represents the general interaction of a two-level atom with an electric 



field in the semiclassical theory, takes the form 

s\{t) = -LUoS2{t) + -[/x, •£(*)] .S3 (i), (1.17a) 

S2{t) ^ LJoSiit) + j:[f^r- E(t)] S3(i), (1.17b) 

Mt) = -j:[f^r ■ E(t)] S2(t) -J:[^^^■ E(t)] Slit). (1.17c) 

It is well known that the dynamical evolution of an A^-level quantum system can be described 
by the rotations of a coherent vector [RRS54, FVH57] which is constrained by the existence of high- 
order constants of motion [ElgSO, HE81, HE82]. For the two-dimensional Hilbert space, some of 
these conservation laws come directly from the properties of the Pauli matrix operators, reflecting 
the intrinsic symmetry of the vector space. For example ct^ = cr| = cr| = /. Another important 
constant of motion is derived from the fact that 

— ](7i{tf = 0, so that al{t)=al{0)=I. (1.18) 

Extending all these ideas to the system of equations (1.17), we may obtain their associated conser- 
vation law, which is written as 

slit) + slit) + slit) = I. (1.19) 

This expression means that the probability of the system is conserved over the time, or equivalently, 
the state of the atom remains normalized in time. 

For electric dipole Am — transitions, we can adjust the arbitrary phases of the connected levels 

so that /jLi vanishes. By denoting 

2 - 2riE 

-tir-E = ^=n, (1.20) 

we can write the semiclassical atomic equations (1.17) in the simpler form: 

Slit) = -iOoS2it), (1.21a) 

S2(t) - Luosiit) + nit)s3{t), (1.21b) 

hit) = -0(t)s2(t). (1.21c) 

Because these equations are the electric-dipole analogues of equations of a spin-1/2 magnetic dipole 
undergoing precession in a magnetic held [Blo46], the vector s(t) is called the electric-dipole "pseu- 
dospin." 



The physical meaning of the expectation values si(t), S2{t), and S3(i) can be interpreted as 
follows. From Eqs. (1.15) and (1.9), it is clear that 83(1) represents the internal energy of the atom 
in ihoJo units; and from Eqs. (1.15) and (1.10), we see that si(i) and S2{t) are both manifestations 
of the atom's dipole moment operator [FVH57, AE75]. 

The pseudospin Eqs. (1.21) can be rewritten as if they were the equations for the precession of 
a solid body upon which a known torque N^ is acting. The superscript F stands for the coordinate 
system of fixed unit vectors |1), |2), and |3). Thus the set of three equations (1.21) can be expressed 
as the single equation: 

s(t) = N^(i) X s(i), (1.22) 

where the vector s has components si, S2, S3, and the torque vector N^ has components 

Nf (i) = -n{t), (1.23a) 

Nf (t) = 0, (1.23b) 

Nf (i) = cjo. (1.23c) 

We note that the pseudospin precession is originated by the first and third components of the torque 
vector. 

To simplify the mathematics, we define a coordinate reference frame which rotates at the same 
frequency lu of the field. In this way we reduce the number of rapidly oscillating variables of the 
system and consider only those which change slowly with time. The torque vector is then rewritten 
as the sum of three torques, one N° along the |3) direction, and two much smaller torques that lie 
completely in the |1)-|2) plane: 

N^ =N+(t)+N-(t)+N°, (1.24) 

where 

N" = (0, 0, cjo), (1.25a) 

N" = (-0 coswi, -fi sincjt, 0), (1.25b) 

N+ = (-0 coswt, +n sinujt, O), (1.25c) 

and 

E{t)^£{t)[e"^* + c.c.]. (1.26) 
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As we see, N"*" rotates counterclockwise as t increases, while N~ rotates clockwise. In a coordinate 
system following s and moving to the right at angular velocity w, the vector N"'' remains constant, 
and N^ is counter-rotating at angular velocity 2uj. In such a coordinate frame the effect of the 
torque N+ on a spin is steady and cumulative over long times. On the other hand, the effect of the 
torque N^ reverses itself 10^^-10^^ times/s, and is almost completely ineffective [Shi63]. 

The rotating-wave approximation (RWA) consists of ignoring N^ for this reason, and writing 
the pseudospin equations using N+ and N^ in place of N^ [BS40, EWG76]. It then follows that 



si ~ —ujqS2 — r2s3 sinwi, 

•S2 = i^osi + r2s3 coscjt, 

S3 — — 5l[s2 coscjt — Si sincjt] . 



(1.27a) 
(1.27b) 
(1.27c) 



By introducing an appropriate rotation matrix for the vector s and defining a nearly stationary 
vector p in the rotating frame with components m, v, and w, it is possible to determine what the 
observer sees in the rotating frame. Thus 



w 



cos Lot sinwi 

— sinwi coswi 

1 



(1.28) 



The equations of motion obeyed by the components of the pseudospin p in the rotating frame 



are 



ii = — (wo — ^)'>^i 

i) — +{ujQ — ^)u + ^w, 

w — —flv. 



(1.29a) 
(1.29b) 
(1.29c) 



which are the same as the single vector equation 



d^'' = ^"^' 



(1.30) 



if the rotating frame torque vector N has the components 



N= {-n, 0, LUo-oj). 



(1.31) 



1.3 TT-pulses 

By defining a "dimensionless" quantity 6{t) as 



t 



d{t) = / n{t')dt', (1.32) 

equations (1.29) can be integrated to give 

u{t;0) = Wo, (1.33a) 

v{t; 0) = Wo sin 0{t) + vq cos 0{t), (1.33b) 

w{t;0) — —vosin9{t) + wocos9{t), (1.33c) 

where uq = -^(0; 0), and so on. The second zero in the labels w(0; 0) and w{0; 0) makes reference to 
the detuning frequency A — loq ^ io. 

In the special case when the applied field envelope has a steady state value between ii and ^2, 
Eq. (1.32) can be integrated to give 

e^no{h-h), (1.34) 

where flo is called the Rabi frequency on resonance. 

The Rabi frequency gives the rate at which transitions are coherently induced between the two 
atomic levels. If the atom is initially in the ground state (wq — ~1, vq — 0) then after a time dt 
such that fl()St — tt, Eq. (1.33c) shows that w = +1, and the atom is in its upper state. In other 
words, a "tt pulse" of electromagnetic radiation inverts the atom population from the ground state 
to the excited state. Now, in the spins terminology, the tt pulse turns a spin from alignment to 
anti-alignment with a static magnetic field. 

The quantity Oo(t2 — ii) is exactly the area under the curve pulse amplitude-time, and represents 
the well-known "area theorem" written as 

A{t)^ f n{t')dt' = 9{t). (1.35) 

Resonant pulses with areas tt, 2tt, 37r, and so on, invert the atomic population one, two, three, and 
so on, times. 

We can also note that solutions (1.33) are the result of a rotation. If the rotating frame Eqs. (1.29) 



are written also as a single vector precession equation 



dt 



p = N X p 



(1.36) 



where 



and 



p = (u, V, w), 



N = {-n,0,A), 



(1.37) 



(1.38) 



then Eq. (1.36) can be represented as 



d_ 
di 






-A 







u 


A 





n 




V 





-n 







w 



(1.39) 



1.4 The three- level atom 

Consider the interaction of a two-mode radiation field with the three-level atom shown schematically 
in Fig. 1.2. Let \a), \c), and \b) represent the initial, excited, and final states of the atom in a 
A configuration. They are eigenstates of the unperturbed part of the Hamiltonian Hq with the 
eigenvalues fiuia, ^c, and fiwf,, respectively {uja < tot < Wc)- 




Figure 1.2: The three-level atom in a A-configuration. 

We are only interested in electric dipole transitions between these three levels, so we can concern 
ourselves exclusively with the three-dimensional Hilbert space spanned by these eigenstates. Then, 
the state vector of the system (in the Schodingcr picture) can be written as 



|V(t)) - C„(t)e-^"»*|a> + Cfc(i)e-^"^*|6) + C,(i)e-*"=*|c), 



(1.40) 
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where Ca, Cb, and Cc arc the slowly varying amplitudes of finding the atom in states \a), |6), and 
|c), respectively. The corresponding time-dependent Schrodinger equation is 

ihj^m)) = Hit)m)), (1.41) 

with 

H{t) = Ho + H,{t), (1.42) 

where Hq and Hi represent the unperturbed and interaction parts of the Hamiltonian, respectively. 
By using the completeness relation 

\a){a\ + \b){b\ + \c){c\^I, (1.43) 

we can write Hq as follows 

Ho = IHoI 

^ huja\a){a\ + hujb\b){b\ + fiwc|c)(c|, (1.44) 

where we have used HqIj) — fiwjjj), with j — a, b, c. Similarly, the part of the Hamiltonian Hi that 
represents the interaction of the atom with the radiation field is described by 

/7i(t)--/x-E(ro,i), (1.45) 

where fi is the atom's dipole moment operator, and E(ro,i) is the two-mode electric field operator. 
In the dipole approximation, the electric field is evaluated at the dipole position Tq ~ 0, and the 
operator can be written in the form 



E(0,t) = ^Re[efe£:fc(t)e'"'=*], (1.46) 

fc=i 

where e^ and £k (t) represent the normalized polarization vector and the amplitude of the electric 
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field, respectively. We can also write Hi in Dirac notation by using Eq. (1.43) 

ffi(t) = -/[/x-E(0,t)]/ 



i,j k—1 



£fe(t)e*"'=*+£:^(i)e-*"'=* 



|J>01- (1-47) 



Assuming that the electric field is linearly polarized along the dipole moment direction, we have 

{i\fi ■ iklj) ^ {i\^i\j) = ^l^J, (1.48) 

and Eq. (1.47) becomes 

Hi{t) = --^Y.^'^A£k{t)e^"'' + £me-'^'']\^){J\, j^a,b,c. (1.49) 

i,j k—1 

In the general case the dipole matrix elements are complex numbers that might be written simply 
as (see 1.6 for an alternative notation) 

My = |Mii|e'"- (1-50) 

In this problem we consider only the relevant dipole transitions \a) -H- \c) and \b) -H- |c), coupled 
by the electric fields Ei and E2, respectively. The transition \a) ^ \b) is forbidden by the dipole 
selection rules. Thus, the only matrix elements of the electric dipole moment that survive are: iiaa 
fibc, with their respective complex conjugates. Plugging Eq. (1.50) into Eq. (1.49), and considering 
only the allowed dipole transitions, we get 



+ |Mfcc|f2e^("^*+«|6)(c| + \ficb\£2e'^^''-^'>\c){b\+U.c.] (1.51) 



We next derive the equations of motion for the probability amplitudes Ca, C'b, and C'c- By 
introducing Eq. (1.40) into Eq. (1.41) and multiplying the resulting equation from the left by (a|, 
we find that 

Ca^^ l]i(t)e^"e'('^i-'^-)*a, (1.52) 



12 



where the Rabi frequency i7i(i) is defined as 



f7i(i) = ^^#^, (1.53) 



and ujca — ^c ^ ^a is the atomic transition frequency. In deriving Eq. (1.52), we have ignored 
counter-rotating terms proportional to exp[±i(a;i +a;ca)] on the right-hand side in the rotating-wave 
approximation (RWA). Similarly, by multiplying instead by (6| and (c|, we find 

C'b^- f72(i)e''^e*^"^~"-'')*Cc, (1.54) 

Cc = ^ f7t(t)e-*"e-*("i-'^<='''*Ca + ^ n';it)e-"^e-'^'^^-'^-''^*Cb. (1.55) 



with the Rabi frequency 



h 



and the atomic transition frequency ujcb — tOc — ^^b- Introducing the detuning factors 

Ai = oji -ojca, (1.57a) 

A2 = W2 - Web - Ai, (1.57b) 

the coupled Eqs. (1.52), (1.55), and (1.54) then reduce to the set 



a = -r!i(t)e'"e'^i*Ce, (1.58a) 

Cb^^ n2{t)e'^e'^^'+^^^'Cc, (1.58b) 

Cc^\ ni{t)e-'"e-'^''Ca + ^ n;{t)e-'^e-'^^'+'^'^'Cb. (1.58c) 

By making a change of variables and choosing values for the phases (see Appendix 4.3.3), it is 
possible to eliminate the exponentials from Eqs. (1.58). It then follows that 

iCa{t) - AiCait) + ^^ CAt), (1.59a) 

iCbit) = (Ai + A2)Cb(t) + ^ Ce(i), (1.59b) 

^CAt) = ^ Ca{t) + ^ Cb{t). (1.59c) 



These are the equations of motion for the probability amplitudes of the three-level A-system shown 
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in Fig. 1.2. The Schrodinger equation for these ainpHtudes in the rotating-wave approximation reads: 



zh-Cit) = H{t)Cit), 
at 



(1.60) 



where 



H{t) 



2Ai fliit) 

2(Ai+A2) n2{t) 

ni{t) n*{t) 



(1.61) 



andC(t)-[C„(t),C6(t),Cc(t)] 

For the case of two-photon resonance (Ai = A2 = 0) and real matrix elements of the dielectric 
dipole moment, Eq. (1.61) takes the simple form 



Hit) ^ - 



Qiit) 

02(t) 

fllit) 1^2 (t) 



(1.62) 



It is easy to verify (see Appendix 4.3.3) that the following linear combination of bare states \a), |c), 
and \b) are eigenstates of the instantaneous RWA Hamiltonian 



\W+) = ^[sin$(t)|a) +cos$(t)|6) + |c)], 
v2 

\W") = cos$(i)|a) - sin$(i)|6), 
\W-) - i=[sin<&(t)|a) +cos$(i)|6) - |c)]. 



where the (time- varying) mixing angle $ is defined by the relationship 



tan$(i) 



^i(t) 

n^it)' 



(1.63a) 

(1.63b) 
(1.63c) 



(1.64) 



When combined with the related photon numbers in the two radiation fields, the eigenstates given 
by Eqs. (1.63) are called the "dressed states" of the matter-field system. Although we do not keep 
track of the photon numbers, we use this name here as well. The (time-dependent) dressed-state 

eigenvalues are 

1 r- „ 1 r- 

(1.65) 



LU+ = +-^Ql + ni, Lu" = o, uj- = — ^Qi + a 
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1.5 The adiabatic following 



An alternative method for population transfer between two states is based on sweeping the pulse 
frequency through a resonance. If the sweep is sufficiently slow, then it is possible to produce 
complete population transfer between the two states that are connected by the resonance [SBK+92]. 
The adiabatic process can be characterized by a steady state process [Mor64]. That is, the rates 
of change of the varying components of the incident laser fields are assumed to be small enough 
that a quasi steady state is maintained throughout the process [Hio83]. These processes have the 
advantage of being insensitive to pulse area, pulse shape, and to the precise location of the resonance. 
Then they are useful for producing population transfer in an ensemble of atoms that have different 
Doppler shifts and different dipole moments. 

The condition for exact adiabatic following of a system from state \a) to state — 16), without 
populating an excited state (which normally undergoes spontaneous emission), can be seen from 
Fig. 1.3. 




-16) 



Figure 1.3: Graphic representation of the Hilbert space for the three-level system in the basis 
{|a), \b), |c)}, and in the basis of the dressed states {IM^"*"), \W'^), and |T4^~)}. The dressed vectors 
are free to rotate in time following the evolution of the system. 

In stimulated Raman scattering processes (STIRAP) the pulses are applied in counterintuitive 
way, that is, the fl2it) pulse precedes the pulse ^i{t), though they partially overlap. In other words 






0, lini 



n,it) 



t^+oo fl2{t) 



(1.66) 



15 



Assuming that 



..a„*m.M), ,i.„, 



then we have 

$(-oo) = and $(+oo) = 7r/2. (1.68) 

Hence the adiabatic state \W'^) coincides with the state \a) before the excitation and with state —\b) 
after it, so that initiaUy only state \W'^) among the adiabatic states is populated. If the excitation 
is adiabatic, then the system will remain in this adiabatic state all time and the population will 
eventually be completely transferred to state \b). Moreover, no appreciable population will reside in 
the intermediate state |c) at any time which makes the transfer efficiency insensitive to decay from 
this state to other states [VS96]. Therefore in the adiabatic process, the evolution of the wave vector 
|^(t)) follows closely the evolution of the "dressed state" |VF''), which goes from a direction parallel 
to I a) to a direction antiparallel to \b). 

We now consider the conditions under which the system evolves adiabatically. Nonadiabatic 
coupling between the eigenstates is small when the rate of change of the mixing angle $(t), is small 
compared with the separation of the corresponding eigenvalues [Mes99] . 

\{W+\W°)\<t:\uj^ ~uj°\. (1.69) 

For no detuning, this separation is given by 

\io^-io''\ = ^^ni + ni^n,s- (i.70) 

In addition, it is easy to find that |(VK+|VF°)| — — 'I'sinG, and therefore the adiabaticity constraint, 
with sin = 1, reads 

|$|«rJeff. (1.71) 

By differentiating the above expression with respect to time, we may get 



nj + ni 



« n,s. (1.72) 



Finally, for a given counterintuitive sequence of pulses ili and fl2, separated by some time 
delay At/T (with T being the pulse length or interaction time), the adiabatic theorem leads to the 
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condition 

f^offT>l. (1.73) 
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Chapter 2 



Laser-induced population transfer 
in multilevel systems 



Cavity quantum electrodynamics (QED) is the part of physics that studies the interaction of single 
atoms and photons coupled to an electromagnetic resonator. Many interesting effects have been 
observed during the last twenty years. They include the alteration of the atomic radiative rates 
by the presence of a cavity around an atom, shifts in the atomic energy level due to the coupling 
with the cavity, manipulation of photons by using the interaction with individual atoms crossing the 
cavity, creation of non-classical field states, among others [OBRW96]. 

In the recent few years, and due to the advancements in the technology of high-Q cavities and 
in atomic beam manipulation, many proposals in the area of quantum information and quantum 
computation have been made. They include applications to particle teleportation, quantum cryp- 
tography, spectroscopy and conditional dynamics. 

In this chapter we study numerically the transfer of population in a four-level system by using 
TT-pulse methods and adiabatic passage schemes. We present some examples of the resonance- like 
features in the failure probability p. We show that the features appear for Gaussian, hyperbolic 
secant, and Lorcntzian pulse profiles. 

2.1 The atom-cavity system 

Consider a three-level atom consisting of two ground states \gi) and jgj), and an excited state 
|e), interacting with a coherent driving field 0(f) of frequency lo^ and a cavity- mode field g{t) of 
frequency w, as indicated in Fig. 2.1. 
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2g{t)^/^^+l 



Figure 2.1: A Three- level atoni in a cavity. 
The state vector describing the three-level atom can be written in the form 

\^{t))=Y.[Cg,At)\gi,n) + CeAt)\e,n)^C,^At)\g2.n)], (2.1) 

n 

where Cg^.n, Ce^m a-nd Cg^.n are the probability amplitudes of finding the atom in states |g]^,n), 
|e, n), and |g2, n), respectively. The ket |gi, n) = |g]^) (g) |n) represents the state in which the atom is 
in the ground state \g^) and the cavity field has n photons. Similar descriptions exist for the other 
two states. 

The time evolution of the system is described by the Schrodinger equation 

ih\<i>{t))^H{t)\^{t)) (2.2) 

with the Hamiltonian operator Hit) given by 

H{t) — tiioa^a + u;eg|e)(e| — ihg{t)(^\e){g2\a — H.c.) 

+ ihn{t){\e) (gi |e-*"^* - H.c.) . (2.3) 

Here a represents the annihilator operator for the cavity mode. The time dependence of ^{t) and 
g{t) may be provided simply by the motion of the atom across the laser- and cavity-field profiles. 

2.1.1 Dark state 

It is interesting to observe that the interaction part of the Hamiltonian (2.3) can only cause transi- 
tions between states within the family 

{|gi,n), |e,n), |g2,n+l)}, (2.4) 
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Therefore, in the rotating-wave approximation, the adiabatic energy eigenvalues of the Hamiltonian 
associated with a particular family of states are [PMZK93] 

En = nfvjj, (2-5) 



E^ = n^w + - [a ± v/A2+4g(i)2(n + l) + f7(i)2j , (2.6) 

where we have assumed that uj = lu^, and A = ujeg — cj is the detuning. We are interested in the 
eigenstate corresponding to En = nhv, which is given by 



1^ .^^^ ^ 2g{t)V^rTl\g^,n}+n{t)\g^,n+l} 
Vf^(t)2+4g(t)2(n + l) 

This eigenstate is, at all times, free of any contribution from the excited state |e, n), and is indepen- 
dent of the detuning A. 

2.1.2 Adiabatic following 

If we assume that only level \gi,n) is initially populated, complete population transfer occurs 

if [KGHB89] 

Oft') p-(t\ 

= 0, (2.8) 



n{t) 



g{t) 



= and ^^'^ 






t— >+oo 

where t — > — oo and t — > +oo corresponds to times before and after the interaction with the fields, 
respectively. Consequently, for the pulse sequence in which the pulse n{t) is delayed with respect to 
g(t), the dark state is the appropriate vehicle for transferring population from state \gi,n) to state 
|g2, n + 1) without populating state |e, n). 
If the condition for adiabaticity [Mes99] 

^oT, 2goV^mr>l (2.9) 

is satisfied (with i7o and gp representing peak intensities for the respective fields) then the state 
vector of the system remains very nearly an eigenvector of the time-dependent Hamiltonian at all 
times, that is 

|*(t)) « \En{t)}. (2.10) 
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2.1.3 Master equation 

We can generalize the previous ideas to more complicated and realistic atomic-level structures. For 
example, this time we may consider Zeeman ground state levels, and include in our analysis the 
spontaneous emission F, and the cavity decay k. 

Consider the case of an atomic Jg = 1 ^ Jg = transition in a four-level A atom, as indicated 
in Fig. 2.2. The transition \gi, 1) -^ |e,0) is strongly connected to a 7r-polarized cavity-mode field 

|eO) 




|g-i,0) |g-o,l) la,0) 

Figure 2.2: A four- level atom in a cavity. 

with coupling strength g{t). The transition \g_i,Q) — > |e, 0) is coupled to a coherent cr^-polarized 
laser field with frequency ljl and Rabi frequency Q,{t). The excited state is assumed to be able to 
decay to the three ground states with the same decaying rate F/3. 

In the short period of time compared to the natural decay times, the dynamical evolution of the 
density matrix p{t) of the atomic system is given by the Liouville equation 



-^ = -i[i7offP-H.c.] +F V A„pA\ + Kapa\ 



(2.11) 



where p{t) is the reduced density operator of the system, and 



//eff=(A-ir/2)^|Jeme)(J< 



e'Tlel - i-a'a 



m{t){A^ 



A 



+iy 



ig{t){a^Ao-Ala), 



(2.12) 



is the effective non-Hcrmitian Hamiltonian in the rotating- wave approximation [PMZ+95]. For the 
two-photon resonance problem (A = 0) this Hamiltonian reduces to 



HcS = - i-^|e)(e| - i-a)a 



in{t){A+i - A\^) - zg(t)(aUo - Ala), 



(2.13) 
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where the atomic lowering operators A^ are given by 

^o- = ^ \Jgmg){Jgmg;la\Jeme){Jeme\, (2-14) 

with [JgiUg] lcr|Je?ne) the Clcbsch-Gordan coefficient for the dipolc transition |e) — > |g) with polar- 
ization (T — 0, ±1. Working out these coefficients, we may have 

A, = ^\g_,){e\, A, = -^\g,){el A^, = ^\g,){e\. (2.15) 

Upon substitution of Eqs. (2.13) and (2.15) into Eq. (2.11), we find that 

|j = -i(r|e)(e|p + «atap + H.c.)-^(|g_i)(e|p-|e)(g_i|p + H.c.) 
- --/^(a^|go>(e|p- |e)(go|ap + H.c.) + napa) 
+ ^(|gi)(e|p|e)(gi| + |go)(e|p|e)(go| + |g_i)(e|p|e)(g_i|). (2.16) 

We can sec from Eq. (2.16) that only the following family of eigenstates survive: 

|g_i,0), |e,0), |go,l), |go,0), and |gi,0). (2.17) 

This means that in our density matrix approach, we have up to 25 matrix elements of the form 
Pij = {Ap\i)- However, we can get rid of some of these matrix elements because they are electric 
dipole forbidden or make no physical sense. In this way, the total number of matrix elements can be 
reduced from 25 to only eight, which corresponds to the number of differential equations describing 
the evolution of the system. 

By introducing the following notation for the atomic levels: 

|a)-|g_i), |6) = |go), |c) = |gi); (2.18) 

and calling p^ = Re[(i|p|j)l , we can write the equations of motion for the density matrix elements 
in the form 

r 

PaOaO = -2 VL(t) paOeO + ^ PeOeO , (2.19a) 

r 

PfcOfcO = -r PeOeO + KPfclfcl, (2.19b) 

22 



Pbibi = -upbibi - 2g(t)ptieo, (2.19c) 

r 

PcOcO = -rPeOeO, (2.19d) 

PeOeO = -rpeOeO +2ri(t)paOeO + 2g(i)pble0, (2.19e) 

P 

PaOeO = -^ PaOeO - ^it){peOeO - PaOao) + g{t) PaObl, (2.19f) 

PaObl = --Z PaObl - ^{t) PbleO - g{t) PaOeO, (2.19g) 

PbleO = --^{T + KJPbleO +^it) PblaO - g(t) (peOeO - Pblbl) ■ (2.19h) 



2.2 Numerical results 

In this section we have studied numerically two different methods for transferring the population 
between two atomic levels: the 7r-pulses and the adiabatic passage methods. The calculation of the 
transfer efficiencies was done by numerically integrating the system of equations (2.19) using the 
Rungc-Kutta method of fifth order with an adaptive mesh. By using Gaussian, hyperbolic secant, 
and Lorentzian pulse shapes, we obtained the coherent transfer efficiencies for the four-level system 
in a cavity, including radiative decay from the excited states. For simplicity, we have considered a 
high-Q cavity with no decay (k — 0). Five different parameters were considered in our simulations: 
the Rabi frequency and the width of the pulses, the time delay (for intuitive and counterintuitive 
configurations), and the spontaneous emission of the excited states. 

Because of the big volume of data manipulated during the simulations (necessary to obtain results 
with some degree of accuracy) the computations were performed in a cluster of four heterogeneous 
computers running Parallel Virtual Machine (PVM). Surfaces of the logarithm of the failure proba- 
bility were plotted for some physically reasonable intervals for the Rabi frequencies and widths of the 
pulses. For example, most of the time we studied widths in the interval [0, 10]. We also explored time 
delays in intervals between [0, 50] for Gaussian pulses, and [0, 100] for Sech and Lorentzian pulses. 
The figures obtained for all these parameters were given in "absolute" numbers (no units). We did 
it in this way because we wanted to examine the ratios between quantities with same dimensions 
rather than individual values taken by the parameters themselves. This is why most of the graphs 
were plotted in terms of fl/g, and At/a. During the numerical simulations we also fixed the value 
g to unity for many of them. This allowed us to have a point of comparison between the pulses' 
energies (Rabi frequencies) and widths for the different failure probabilities obtained, and then to 
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conclude which set of parameters achieved the maximum transfer efficiency. 

Our analysis starts with the 7r-pulses method. An intuitive sequence of these pulses can, in 
principle, produce a complete population transfer between two atomic levels [SBK+92]. Here, we 
adjust the laser intensities and pulses' duration so that the time integral of the Rabi frequencies 
(the pulse area) have the value of tt. The first pulse 17, takes the system from the initially populated 
level \g_i 0) to the excited level |e 0); then the cavity- mode field g takes the system to the target 
level |go 0). For this to be possible the excitation has to be coherent. Results for the time delay and 
failure probability are shown in Table 2.1 for different values of the Rabi frequency and the width 
of the pulses. 

Table 2.1: 7r-pulses method. 



Pulse profile F O a g ag At logj^pp 

Gaussian 0.01 2.14 0.29 1.00 0.63 1.26 -2.05 

0.02 2.15 0.29 1.00 0.63 1.16 -1.78 

0.05 2.03 0.31 1.00 0.63 1.01 -1.44 

0.10 2.11 0.30 1.00 0.63 0.89 -1.19 

0.20 2.18 0.29 1.00 0.63 0.75 -0.95 

Sech 0.01 1.84 0.26 1.00 0.50 1.52 -1.94 

0.02 2.01 0.25 1.00 0.50 1.35 -1.69 

0.05 2.16 0.23 1.00 0.50 1.11 -1.37 

0.10 2.26 0.22 1.00 0.50 0.94 -1.13 

0.20 2.43 0.21 1.00 0.50 0.77 -0.91 

Lorentzian 0.01 5.71 0.09 1.00 0.50 2.28 -1.63 

0.02 5.53 0.09 1.00 0.50 1.80 -1.43 

0.05 5.56 0.09 1.00 0.50 1.32 -1.16 

0.10 5.80 0.09 1.00 0.50 1.04 -0.96 

0.20 6.28 0.08 1.00 0.50 0.82 -0.77 

As is well known, this method presents the most efficiency (100%) when the system undergoes no 
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spontaneous emission. However, when the system decays radiatively, the efficiency of the method 
decreases graduahy as the spontaneous emission rate increases. We observe that, when the sponta- 
neous emission is F = O.Of , the Gaussian pulses have an efficiency of 99.1%, while the secant and 
Lorcntzian pulses presented maximum efficiencies of 98.9% and 97.7%, respectively. These values 
were the best obtained for a wide range of the examined pulse amplitudes and pulse widths. Fig. 2.3 
shows that for greater values of the spontaneous emission, the transfer efficiency reduces in a consid- 
erable way. The time delay between the pulses also reduces to shorter times for bigger spontaneous 
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Figure 2.3: 7r-pulses failure probability for four different values of spontaneous emission. The values 
of the Rabi frequencies, widths, and time delay can be found in Table 2.1 

emission values. This is understandable because more overlap between the pulses is required to 
neutralize the decay rate of the excited states. 

It is interesting to note that in this intuitive sequence of pulses (as it is the case for the 7r-pulses 
method) , good transfer efficiencies were achieved when the driving classical field (laser) was stronger 
than the quantized cavity mode, as shown in Fig. 2.4. Even though the ratio of the Gaussian pulses 
(^/g) was bigger than the ratio of the Sech pulses, less interaction time {/St/ (jg) was required for 
the Gaussian pulses to obtain a bigger efficiency. The Lorentzian pulses required the biggest energy 
(Rabi frequency) to achieve even less efficiency than the other two types of pulses, as observed from 
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Figure 2.4: 7r-pulses' profiles for tfie best set of paranieters when F = 0.01. 
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the ratio of the pulses shown in 2.4. 

The TT-pulses method is a good and efficient technique to transfer the population in atomic 
systems, but we will see that the adiabatic passage scheme overtakes considerably the 7r-pulses 
performance for some very interesting and particular situations. 

The next method we studied numerically was the adiabatic passage scheme. Table 2.2 shows 
the values of time delay and failure probability we obtained for different values of the spontaneous 
emission. 

Table 2.2: Adiabatic passage scheme 



Pulse profile F il a g (jg At logioP 

Gaussian 0.00 2.00 1.00 2.00 1.00 1.31 -4.88 

4.00 1.00 19.20 1.00 1.90 -4.53 

6.00 1.00 5.70 1.00 1.50 -6.83 

6.70 1.50 2.00 1.00 2.72 -5.67 

0.10 3.39 3.23 1.00 2.45 5.85 -1.99 

2.75 3.09 1.00 2.48 5.29 -2.00 

0.20 3.30 3.30 1.00 2.50 5.90 -1.73 

2.40 3.30 1.00 3.00 5.38 -1.80 

2.30 4.20 1.00 4.48 6.62 -1.96 

2.10 4.60 1.00 5.00 7.09 -2.01 

Sech 0.00 2.00 1.00 2.00 1.00 0.80 -7.79 

0.10 2.60 1.40 1.00 1.20 2.70 -1.71 

5.00 1.50 1.00 2.00 3.80 -1.90 

4.50 1.70 1.00 2.40 4.10 -1.89 

4.20 1.80 1.00 2.60 4.30 -1.86 

6.30 3.40 1.00 4.01 12.01 -2.00 
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Table 2.2: Adiabatic passage scheme 



0.20 14.70 5.00 1.00 7.00 21.40 -2.06 

Lorentzian 0.00 2.00 1.00 2.00 1.00 0.32 -4.46 

0.10 9.20 0.60 1.00 2.00 2.29 -1.05 

Here, it is important to em.phasize that the figures shown in the table do not necessarily correspond 
to the best possible values of the transfer efficiency achieved by the system. It happens that the 
efficiency can be optimized as much as we want by increasing the values of the Rabi frequencies, or 
equivalently, the widths of the pulses. Thus, for bigger Rabi frequencies, bigger efficiencies. This 
is something that has no experimental worth. What we want to examine here are those values of 
the Rabi frequency and time delay that can be reproduced in a laboratory, and for which very good 
transfer efficiency values can be obtained. So we have decided to consider a reasonable value of 
logioP = —2.00 for our numerical exploration. 

We analyze first the case of no spontaneous emission, when the system is driven by Gaussian 
pulses. Fig. 2.5 shows that for the particular set of parameters {fi = 2.00, a = 1.00, g = 2.00, ag = 
1.00, At ~ 1.31} a tremendous efficiency of 99.9% or higher is obtained. We observe that, while the 
failure probability remains almost the same over a wide region of the pulses' amplitudes, there is 
one particular sector for which a sudden fall in the failure probability occurs. A very sharp, deep, 
and unexpected well just appears on the probability surface. Intriguingly, an even more sharper 
and deeper well emerges when the driving pulses are hyperbolic secants, and the set of parameter is 
{n = 2.00, a == 1.00, g == 2.00, ag = 1.00, Ai = 0.80}, as shown in Fig. 2.6. For this case we observe 
that the region nearby to the sharp dip is deformed a little bit forming a valley of good values for 
the transfer efficiency. Fig. 2.7 shows another efficient set of parameters, this time for Lorentzian 
pulses. 

We observe from Fig. 2.8 that, with a shorter interaction time {At/dg) than the Gaussian pulses 
(but longer than Lorentzian pulses), the Sech pulses achieved the best transfer efficiency. Here, all 
the profiles have the same ratio for the pulses {il/g= 1.0). 

When we consider spontaneous emission, the situation changes a little bit. We are looking for 
good efficiencies of the order of log^gp = —2.00. Figs. 2.9 and 2.10 show the failure probability 
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Figure 2.5: Gaussian surface and contour for F = 0.00. Here fl = 2.00, cr = 1.00, g = 2.00, a-g = 1.00, 
At = 1.31, and log^oP = -4.88. 



29 



Io9io (P) 




Ai/a-f 



(a) Surface 



logio (p) 
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Figure 2.6: Hyperbolic secant surface and contour for F = 0.00. Here il — 2.00, cr == 1.00, g~ 2.00, 
cTg = 1.00, At = 0.80, and log^op = -7.79. 
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Figure 2.7: Lorentzian secant surface and contour for F = 0.00. Here fl = 2.00, a = 1.00, g : 
cTg = 1.00, At == 0.32, and logioP = -4.46. 
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Figure 2.8: (a) Cross-sections of log^gp as a function of At/tTg for F = 0. (b)-(d) Pulses' profiles. 
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for two different values of the spontaneous emission, when Gaussian and secant pulses are used for 
driving the system. In this case, we focus our attention in the Rabi frequencies, the pulses' widths 
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Figure 2.9: Adiabatic passage for F = 0.10 



and the time delay. We immediately observe from the figures that a bigger energy and a longer time 
delay is required by the Sech pulses to obtain the same efficiency than that achieved by Gaussian 
pulses. We also see that a considerable overlapping of the pulses is necessary to obtain good transfer 
efficiencies. Gaussian pulses seem to be more effective this time, as the comparison between the 
ratio of the pulses and the interaction time confirms. 

Unfortunately, it was very difficult to obtain results for the Lorentzian pulses for values of the 
spontaneous emission different from zero. These pulses are very wide, and their numerical integration 
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Figure 2.10: Adiabatic passage for F — 0.20 
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is a very complicated task. However, it is possible to see from Table 2.2, that for bigger values of 
the spontaneous emission, bigger values of the Rabi frequency are needed. 
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Chapter 3 



A simple analytical model 



In this chapter we introduce a simple analytical model that explains qualitatively the origin of the 
nonadiabatic, resonance-like features observed in adiabatic passage methods. We also present a 
mathematical description of the three different pulse shapes used in our numerical simulations of 
chapters two and three, including the nonadiabatic coupling terms. By using the analytical model, 
we study the dependence of the failure probability p on the system parameters fl T, showing the 
exact solution of a simple problem. Although the analytical results here reported have been obtained 
by others, the method is somewhat different from those in the literature. 

3.1 Population transfer 

Consider a three-level A system driven by a two-mode classical coherent field, as shown in Fig. 3.1. 




Figure 3.1: Three- level A system driven by two lasers. The Rabi frequencies for the pump and Stokes 
laser are fli and Q2i respectively. 

We want to study the problem of coherent population transfer from the initially populated state \a) 
to the final state \b) by using a pump laser coupling states \a) and |c), and a Stokes laser coupling 
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states |c) and \b) (see [KGHB89]). 

Under two-photon resonance condition and by using the rotating-wave approximation, the evo- 
lution of the system is described by 



'^a — 7: L'c 



n2it) 



a 



(3.1a) 
(3.1b) 
(3.1c) 



where Ca, Cb, and Cc are the probabihty amphtudes of finding the system in states \a), \b), and |c), 
respectively. 

To obtain approximate solutions to this system, we may proceed as follows. Let us first consider 
the normalized state vector |^(i)) in spherical coordinates, as indicated in Fig. 3.2. The state vector 




\ l*(i)) 



a 



\b) 



Figure 3.2: Graphic representation of the Hilbert space for the three-level system in the eigenbasis 
{\a), \b), \c)}. The state vector |^(i)) is represented by its three Cartesian components Ca, Cb, and 
Cc, corresponding to the probability amplitudes. 

components Ca, Cb, and Cc (probability amplitudes) can be written in the form 



Ca = sin^cost/), 
Cb = sin 9 sin (p, 
Cc = cos 9, 



(3.2a) 
(3.2b) 
(3.2c) 
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with time derivatives given by 

Ca = ^ cos 6' cos (^ — (/) sin sin (/), (3.3a) 

6*6 = ^cos6'sin(/) + 0sin6'cos(/), (3.3b) 

C'c = -^sin6'. (3.3c) 

Now, plugging Eqs. (3.2) and (3.3) into Eqs. (3.1), and by using the mixing angle 

.a„*m = Ml, (3.4) 

it is simple to verify that Eqs. (3.1) reduce to a system of two coupled nonlinear differential equations 
given by 

tan|"6l-7r/2] , . . , .i , ■. 

cj){t) ^ ^— — '—^n{t) cos [0 + $(t)] , (3.5a) 

^W=-^sin[0+$(i)]. (3.5b) 



Here ^{t) = ^ nl{t) + nl{t) . 

Taking a look at the asymptotic behavior of Eqs. (3.1), we see that at very early times (i — > — oo) 
when the state vector is parallel to the initially populated state |a), the probability amplitude 
coefficients take the values 

Ca(— oo) = 1 = sin0(— oo) cos (/)(— oo), 

C6(-oo) = O = sin6l(-oo)sin0(-oo), (3.6) 

Cc(— oo) = = cos6'(— oo). 

This implies that before the interaction with the lasers, the angular coordinates of the state vector 
are 

6'(-oo) = tt/2 and (/-(-oo) == 0. (3.7) 

Similarly, for very late times (t — ?■ +oo) when the state vector is antiparallel to the final state \b) 
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(see Sec. 1.5), the probability amplitudes are 

Ca (oo) = = sin0(oo) cos(/)(oo), 

Cb(oo) = — 1 = sin0(oo) sin(/)(oo), (3-8) 

Cc(oo) = = cos(^(oo). 

From these equations, we may conclude that after the interaction with the lasers 

9{oo) = tt/2 and (j){oo) = -7r/2. (3.9) 

Since we are interested in transferring population from the initial state \a) to the final state \b) 
without populating the leaking excited state \c) (which we assume undergoes radiative decay), then 
the state vector |^(t)) should evolve adiabatically by closely following the "dressed state" |VF°) of 
the Hamiltonian of system (3.1). Thus the angular coordinates of the state vector can be written as 

9ad{t) = tt/2 and ^ad{t) = -$(t), (3.10) 

where the label "ad" accounts for an adiabatic following of these angles to the corresponding angular 
coordinates of the dressed state. 

If we assume small deviations from conditions (3.10), the system (3.5) can be linearized. By 
plugging the equations 

6 = ead + Se = T:/2 + Se, (3.11a) 

(I) = ^ad + S(l) = -Ht) + S^, (3.11b) 

and their time derivatives 

9 = 59, (3.12a) 

<j) = -$(t) + 6(1), (3.12b) 

into Eqs. (3.5), we obtain 

d(j> = ^69 + $(t), (3.13a) 

S9 = -^S(t). (3.13b) 
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To integrate Eqs. (3.13), we need to make a change of variables by defining a "dimensionless 
phase" r in the form 

(3.14) 



T{t) = I — ^ dt' 



oo 



This phase wih play the role of a rescaled time in the following. Then, by using the chain rule 



d _ dr d _ n{t) d 
dt dt dr 2 dr ' 



(3.15) 



we can rewrite Eqs. (3.13) in the form 



dr dr 

-^6e[T] = -S(t)[T]. 
dr 



(3.16a) 
(3.16b) 



Now, by differentiating the second equation and plugging it into the first, we obtain the equation of 
an undamped driven harmonic oscillator 



clT'^ dr 



(3.17) 



The general solution to this equation is very well known and given by (see Appendix 4.3.3) 



r] =^sinfT] +Bcos[r] - / sin[ t-tM , , dr' . 

Jo dr' 



(3.18) 



The lower limit of the integral results when the time is — oo giving 



T — oo = 



Vt{t) 



dt^O. 



(3.19) 



In addition, we have that the difference between phases can be written as 



r-r^l'^„r 



(3.20) 



Now, returning to the original time scale in Eq. (3.18), we get 



Se{t) = Asin 



^dt' 



n{t" 



B cos 



2ffi.,' 



■dt" 



''^^dt'- 
dt' ' 



(3.21) 
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and by inserting Eq. (3.21) into (3.13b), we can solve for 5(f) obtaining 



S(j){t) = -A cos 



' ^dt' 



i?sin 



' m,, 



cos 



».." 



dt' 



The initial conditions for the linearized system are 



(3.22) 



(56l(-oo) = 0, 

6e{-oo) = -rj(-oo)$(-oo)/2, 
S(f>{-oo) = $(-oo), 
S(f>{-oo) = $(-oo). 



(3.23a) 

(3.23b) 
(3.23c) 
(3.23d) 



And finally, by evaluating the constants of integration A and B, we may obtain 



64>{t) = $(-cx))cos 
6e{t) = -$(-oo)sin 



2 



— oo 

t 



t' 2 



* n{t" 



■ dt'' 



^^^dt' 
dt' ' 

dt' 



(3.24a) 
(3.24b) 



The two terms in each of the equations (3.24) represent different things. The first term represents 
a coherent mixing of the adiabatic states that occurs when the initial state of the system is not one 
of the eigenstates of the interaction Hamiltonian. The second term, on the other hand, represents 
corrections to the adiabatic approximation. 

3.2 Asymptotic behavior of the pulses 

In the present work, we used three different kinds of pulse shapes for the numerical computations 



Gaussian pulses: 0(i) = ilo exp(— i^/2cr^). 

Hyperbolic Secant pulses: 0(t) = JIq sech(t/(T), 
Lorentzian pulses: 0(t) = Oq cP' lit? + cP'^. 



(3.25a) 
(3.25b) 
(3.25c) 



Here JIq is the Rabi amplitude of the pulse, and cr is a scale parameter (width) describing the extent 
of the pulse. The different shapes are illustrated in Fig. 3.3. 

It is important to consider the asymptotic behavior of functions (3.25) because this determines if 
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Gaussian 

Sech 

Lorentzian 



Figure 3.3: Profiles of tire ttiree pulses used in this work. 

a sequence of pulses satisfies the adiabatic conditions, and then the system can evolve adiabatically. 
In addition, the nonadiabatic coupling matrix element 



Ht) 






can be used as a "local" adiabaticity criterion for a given shape of the laser pulses. 

3.2.1 Gaussian pulses 

Consider two Gaussian pulses of the form 



(3.26) 



fliit) = nwexp[-{t - Atf /2af], 



^2(t) = ^20 exp(-tV2cr, 



(3.27a) 
(3.27b) 



These two pulses are in a counterintuitive configuration, where the pulse 2 precedes the pulse 1, and 
their separation is given by the time delay At. 
The ratio of the pulses can be written as 



ni{t) _ ill 



■ exp 



(At)^ 



f^2 ~ cn ,2 



t^ + ^rt 



2ai 2(^1^2)2 a 



(3.28) 



We observe that at very early or very late times, the behavior of this ratio depends on the pulses' 
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widths as follows 



cFi > (72 '■ r2i/i72 — > +00 when t — > ±00, 

when t — > —00, 

-00 when t — > +00, 
(Ti < CT2 : r2i/r22 — > when t — > ±00. 



(Ji = 0-2 : ^i/fl2 



(3.29a) 
(3.29b) 

(3.29c) 



Equation (3.29b) shows that for pulses with equal widths (cti —(72) the adiabatic conditions for 
the system are satisfied (that is, the external field and the coupling field go to zero when the time 
goes to —00 and +00, respectively) and the adiabatic theorem supports the population inversion 
between the initial and the target levels. On the other hand, for unequal widths (tri 7^ (T2) the 
adiabatic conditions are not fully satisfied by the system, and the adiabatic theorem alone cannot 
validate the transferring process. In this case we need additional arguments to help us to understand 
why, in some situations like these, large transfer probabilities are still achieved. Fig. 3.4 shows the 
asymptotic response of the pulses' ratio. 



20 

15 
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,: / 

5- // 






0"! > (T2 
O"! = 0-2 



\ 



V 



-10 -5"' " 5 10 15 time 

Figure 3.4: Behavior of the Gaussian pulses ratio as a function of time, and for different pulse widths. 



The nonadiabatic coupling for these Gaussian pulses is given by 



. ^ [l]iol]2o/(gig2)'] [<r2^t - [al - al)t] 
^ ^ " 0?o exp(tV'Ti) + f^io exp[(t - Mf/al] '""P 



(t - Ai)2 e 



2al 



2ai 



(3.30) 
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3.2.2 Hyperbolic secant pulses: 

Consider two Sech pulses of the form 



ni{t) = riw scch[(i - At)/(Ti] , 

ri2(i) = ^20 SCch(t/(T2)- 



(3.31a) 
(3.31b) 



The ratio of the two pulses is given by 



rji(t) fn 



'10 



exp(i/(72) +exp(-t/cr2) 



^2{t) V^2o/exp[(i- Ai)/cri] +exp[-(t- Ai)/cri]' 
and the asymptotic behavior of the pulses can be expressed as: 

ci > (72 : ni/fl2 — > +00 when t — > ±oo, 

(i7io/i^2o) exp(— At/crj when t — > — oo, 

(i7io/r22o) exp(Ai/crj when t — > +oo, 

(Ti < CT2 : r2i/i72 ^^ when t — > ±oo. 



CTi = 0-2 : rii/r22 



(3.32) 



(3.33a) 
(3.33b) 

(3.33c) 



From these relations we observe that hyperbolic secant pulses behave like Gaussian pulses when 
the widths are unequal (cri 7^ cr2). In contrast, for equal widths (cri = (72), the ratio of the pulses 
reaches constant values different to zero or infinity. In both situations the adiabatic conditions are 
not satisfied by the system. However, we will see in Ch. 2 that Sech pulses still allow us to get very 
good transfer efficiencies. 

For Sech pulses, the nonadiabatic coupling is given by the function 



m 



[Oiori2o/o'i0'2] sech[(t- At)/cri] sech[t/cr2] 
nig sech^ [(t - At)/ai] + fl^g sech^ [t/trs] 

't- At 



a\ tanh 



(^) 



— (72 tanh 



CTi 



(3.34) 
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3.2.3 Lorentzian pulses: 

Finally, let us consider two Lorentzian pulses of the form 



(3.35a) 
(3.35b) 



The ratio of the two pulses is given by 



ni{t) fnioal\ t^ + cFj 



n2{t) 



and the asymptotic behavior expressed as 



no 



'^2 J {t-Aty + aj 



(3.36) 



(Ti > (72 '■ fti/n2 ^^ ^10 o'i/^20 '^2 when t — > ±oo, 

(Ti = (72 : f7i/r22 -^ ^io/f^20 when t — > ±oo, 

CTi < (72 : ili/Jl2 —7- flio ai/n2o 02 when t — > ±00. 



(3.37a) 
(3.37b) 
(3.37c) 



These equations show a very interesting property: whatever the time delay of the pulses is, their 
ratio is the same constant in both directions (±f). Lorentzian pulses behave in a very strange and 
sometimes unpredictable way. Ch. 2 will show some numerical results obtained when Lorentzian 
pulses were used to drive a quantum system. 

The nonadiabatic coupling of the Lorentzian pulses is given by 



m 



2r2iof^20Cr?cr|{t[(t- At)2 + cr2] _ (f - At) (t^ + cr|) } 



(3.38) 



3.3 The failure probability 

To determine the degree of efhciency in the process of transferring atomic population from the initial 
state I a) to the final state |6), we define the "failure probability" as 



p = 1 - |Cb(cx))| 



(3.39) 



This quantity gives us information about how much the process fails to get the population transferred 
to the target state. If we substitute Eqs. (3.2b) and (3.11) into Eq. (3.39), and consider only the 
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lowest order in S9 and 5(j), the probability failure can be written in the form 



P=[{Scf,r + (S9)']^_ 



(3.40) 



Let us first consider the case of two Gaussian pulses driving the three-level system described by 
equations (3.1). If the pulses' widths are considered equal (cri = (72) then the system satisfies the 
adiabatic conditions and the mixing angle takes the values $(—00) = and $(00) = 7r/2. Solutions 
for the small angle deviations (3.24) can now be written as 



se{t) 



n{t") 



f 2 

t 



dt" 



d^{t') 



VL{t") 



dt" 



dt' 
d<^{t') 
dt' 



dt'. 



dt'. 



(3.41a) 
(3.41b) 



Upon substitution of these equations into Eq. (3.40), and by using the relation 



<±i5e\ ={5(t>f + {5ef, 



(3.42) 



we find the following closed-form expression for the failure probability 



P = 



exp I 



^ '■ dt' — dt 
dt 



(3.43) 



It can be readily seen that integral (3.43) is small enough to explain the large transfer probabilities 
observed in adiabatic passage processes. As long as the adiabatic conditions are satisfied {Via 3> 1), 
the exponential term in the integrand can oscillate rapidly over the interval of integration, making 
the integral to take small values. The local adiabaticity criterion also suggests that the nonadiabatic 
coupling d^/dt remains small over that interval, contributing with a small value to the integrand. 
Similarly, we note that the very sharp and profound dips found in the Ch. 2 simulations for some 
particular choices of the pulses' parameters (see Figs. 2.5 and 2.9) correspond to zeros, or near-zeros 
of Eq. (3.43). 

By making the change of variable 



Tit) 



* n(t') 



dt'. 



(3.44) 
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the probability failure (3.43) can be written in the form 



P = 



(3.45) 



-,— ^r 



e ~'——dt 
dt 



This expression can be interpreted as the Fourier transform of a function equal to d^/dr between 
T = and t(oo), and zero elsewhere. This accounts for the oscillatory character of integral (3.43), 
and explains the existence of zeros (or near-zeros) for certain values of the parameters. 

By way of example, let us consider a simple case for which integral (3.43) can be evaluated 
analytically. If the pulses are given by the equations 

l)i(t) = nosm{t/T), 

f7o cos(t/T) if -T7r/2 < t < TTr/2, 
elsewhere. 




we may have that 



$(i) = |, and n{t) = Re- 



plugging these values into Eq. (3.43), we find that 



P 






16 sin2(7rf7or/4) 



nlT^ 



This result reveals a power-law decay of the failure probability in the form l/{flo T)'^, and also shows 
zeros for this function when Q,qT = An. 

Although it is well known that adiabatic passage produces very large transfer efficiencies provided 
the adiabatic conditions are satisfied, the dependence on 17 T of the failure probability is complicated. 
A "perfect adiabatic transfer" is characterized by a probability failure that decreases exponentially 
as a function of OT (this is called the DDP result [DP76]). Recent results have shown that this 
exponential behavior is no longer satisfied for sufficiently large values of il T, giving instead a power- 
law decay of the probability failure with respect to these parameters [LS96, VS96]. In addition to this 
power-law, oscillations of the failure probability have been also observed coming from a nonadiabatic 
component of the total transition amplitude [DII98]. 

Now, let us explore a bit what happens when the conditions for adiabaticity are not fully satisfied, 
for example, the case of hyperbolic secant pulses. As before, we want to transfer population from 
the initial level to the target level by closely following the evolution of the dressed (adiabatic) state 
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of the system. If the pulses' widths are in such a way that the ratio rii/i72 goes to constant values 
at very early or very late times (see Eq. (3.33b)), the adiabatic conditions are obviously not satisfied 
by the system. Then an adiabatic following of the dressed state by the state vector is not possible, 
at least under those initial conditions, because these to vectors are not parallel to each other before 
the interaction with the lasers takes place. In other words, the state vector is parallel to the initially 
populated level, while the dressed vector is somewhere else (see Fig. 1.3 for a representation of these 
vectors in the Hilbcrt space). However, it happens that at some moment during the evolution of 
the system, the dressed state approaches the initial state and the target state close enough to take 
sufficient population from the initial level, evolve parallel to the state vector for a moment, leave 
the population in the target state, and finally continue with its own evolution to some other place 
in the Hilbcrt state. Although the system does not evolve in an adiabatic manner, strictly speaking, 
there is a region in time (not too early, not too late) for which this evolution resembles that of an 
adiabatic process. In this "near-adiabatic" process, large transfer efficiencies are still achieved. 

Fig. 3.5 shows the evolution of the system for two particular set of parameters corresponding 
to Gaussian and Sech pulses. We observe from these figures that the mixing angle is close to zero 
for a short period of time before the interaction, and nearly 7r/2 for a short period of time after 
the interaction. When t -^ — oo, the dressed state |VF°) is parallel to the empty bare state \b) and 
carries no population. However, the dressed "unstable" states \W^) and \W^) (see Sec. 1.5 have 
components along the bare state \a) (which has the initial entire population) and each one carry 
half of the total population. These two dressed states also have components along the excited state 
|c) which decays radiatively The adiabatic frame of reference rotates faster than the frequency 
of the lasers (and perhaps than the decay rate itself) and the coupling classical field J7(i) has no 
populated enough the decaying level. Then the system is not going to lose too much population at 
the beginning of the interaction with the cavity- mode g. By the time in which fl(t) starts to act 
upon the levels, the adiabatic state |VK") is almost parallel to the initial state \a). At this point, this 
adiabatic state will carry mostly the entire population of the initial state (not all of it because the 
other unstable adiabatic states have already contributed with population to the target state) and 
the wave vector will follow adiabatically the evolution of this adiabatic state. 

Our simple analytical model can also explain the behavior of the system for those near-adiabatic 
conditions. Since the goal is to have 

0((^) = -|, and 9{^) = ^, (3.46) 
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Figure 3.5: Spherical angles. For the Gaussian pulses we have fl = 2.75, a = 3.09, g = 1.00, 
(Tg = 2.48, Ai = 5.28, and logj^pP — —2.00. For the sech pulses we have O = 6.30, a ~ 3.40, 
g — 1.00, (Tg = 4.01, At = 12.01, and logj^pp = —2.00. In both cases the spontaneous emission was 

r = o.io 
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we see from equations (3.24) that the condition 



oo 2 



dt =■ riTT, with n odd (3.47) 



must be satisfied. This ensures that the sine in the first term of the S9 equation cancels out. The 
second terms of equations (3.24), which are proportional to $, can be ignored for a moment (we 
will retake these components later on). Therefore Scf) ~ —$(—00). By plugging this result into 
Eq. (3.11b), and considering the simple case in which flio = O20, we may obtain 

^(cx)) ~ -$(00) - $(-00) 
~ -[$(00) + $(-00)] 
~ -[arctan(e^*/'^) + arctan(e-'^*/'')] 
~ -7r/2. (3.48) 

Now, we have to examine what happens with the second terms in (3.24). Since those terms are 
typically small, because of the near-adiabatic condition (that is, they are inversely proportional to 
ri cr), a small deviation from the parameter values that make the first terms in (3.24) vanish exactly 
may be enough to, instead, make the first and second terms (nearly) cancel each other, and this is, 
in fact, what we observe in our numerical calculations. 
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Chapter 4 



Transfer of atomic coherence 



The concepts of atomic coherence and interference has been extensively studied and apphed to 
many areas of atomic physics and quantum optics for the last fifteen years. Applications of these 
ideas include electromagnctically induced transparency (EIT), electromagnetically induced absorp- 
tion (EIA), lasing without inversion (LWI), sensitive spectroscopy in coherent media, among other 
things [FYLOO]. 

Recent developments in quantum information and quantum computation have stimulated interest 
in processes that can prepare entangled states, as well as perform conditional quantum dynamics 
and logic gates. The coherent manipulation of these entangled states are fundamental to realizing 
a quantum computer, and promises a novel atomic spectroscopy with resolution better than the 
standard quantum limit [PGCZ95]. 

In this chapter we present a numerical analysis of the Zeeman coherence transfer between two 
atoms by using adiabatic passage methods. We also present a simplified model for the two-atom 
-|- cavity system that help us to better understand the process of coherence transfer between two 
atoms coupled with a cavity. 

4.1 The two- atom + cavity system 

We consider a system consisting of a single-mode cavity containing two three-level atoms in the A 
configuration, as shown in Fig. 4.1. The atoms are fixed inside the cavity at distances apart much 
larger than the wavelength of the cavity mode and interacting individually with laser beams. The 
transitions \a)j -^ \c)j of both atoms {j — 1,2) are coupled to separate classical coherent driving 
fields with frequency ujl and Rabi frequencies flj. The transitions \b)j — > \c)j are strongly coupled 
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atom 1 atom 2 

Figure 4.1: Transfer of atomic coherence between two A systems. 

to the same quantized cavity- mode field with coupling strength g and frequency uj. The dynamics 
of this system is described by the interaction Hamiltonian: 



^/ = 2 E [f^.We-^"^*|c),,(a| +g|c),,(6|a] +H.c., 



(4.1) 



where a is the annihilation operator for the cavity mode and g gives the coupling strength between 
the atoms and the field mode (vacuum Rabi- frequency). 

4.1.1 Probability amplitude method 

We are interested in transferring the atomic coherence according to 



{A\a)i + B\b)i) |6)2|0)e -^ \b)i {A\a)2 + S|6)2)|0), 



(4.2) 



where A and B are arbitrary coefficients. If we consider the particular case in which the system 
evolves without populating the dark state |&)i|fc)2|0)c (that is, for B = {)) then the evolution of the 
quantum system takes place as follows: 



l«)l|^)2|0)c 
|oi(t) 

|c)l|&)2|0). 



\h)i\b)2\l) 



\b)i\ah\0)c 

02(t) 

l&)l|c)2|0). 



(4.3) 
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The state vector of the system can be constructed as a Hnear superposition of these basis eigenvectors 
in the form 

|*(i)) = Cab{t) \abO)+C,b{t) \cbO)+Cbb{t) \bbl) 

+ Cbcit)\bcO)+Cbait)\baO), (4.4) 



with the probabihty amphtudes given by 

Cab{t) = {abO\'i'{t)), Ccb{t)^{cbO\'^{t)), Cb,(t) = (661|*(t)), 
Cbcit) = (6cO|*(i)), Cbait) = (6aO|*(t)). 



(4.5) 



The corresponding time-dependent Schrodinger equation is 

ih\4fit))=Hi\^{t)), (4.6) 

with 

Hi = ^ \c)n{a\ + f |c)n(&|a + ^ \ch,{a\ + | \ch,{b\a 

+ ^ \a)Mc\ + f a^ \b)n{c\ + ^ \ah,{c\ + f a^ |6>22(c|. (4.7) 

By introducing Eqs. (4.4) and (4.7) into Eq. (4.6), and muhiplying the resuhing equation by 
every bra vector of the basis (4.3), we find that 



(4.8a) 
(4.8b) 
(4.8c) 
(4.8d) 
(4.8e) 



These are the equations of motion for the probabihty amphtudes of the two-atom + cavity system. 
Here we have considered the spontaneous emission of the excited levels by adding phenomenological 
decay terms F to the equations. We have also avoided the complex character resulting in some of 
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the equations by introducing the notation 



C, 



be 



iCi 



be 



and 



C 



be 



iC, 



be- 



(4.9) 



4.1.2 Dark state 

Under two-photon resonance condition, and by using the rotating-wave approximation, the evolution 
of system (4.8) is governed by the non-Hermitian Hamiltonian 
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(4.10) 



This operator has a dark state with eigenvalue zero given by 



\D) = 



n2g\abo) -nin2\bbi) + nig\baO) 

7WF+^2 ' 



(4.11) 



where V,o = \/^'i + ^2- 

We observe that states of the type given by Eq. (4.11) are immune against decay from the excited 
atomic levels, since they have no contribution of such states. Therefore, the dark state \D) becomes 
the right mechanism to transfer, from one atom to the other, the Zeeman coherence of the ground 
state levels by using adiabatic passage [PGCZ95]. 



4.1.3 Conditions for adiabatic following 

If the pulses are applied in a counterintuitive sequence, that is the pulse on atom 2 preceding the 
pulse on atom 1, an adiabatic transfer of the dark state between the atomic levels \abO) and \b aO) 
may be achieved. To see clearly how this happens, let define the mixing angle $ by the relationship 



tan$(f) 



»i(t) 
n2{t)'' 



(4.12) 



and express the dark state (4.11) in the form 



\D)^K 



cos$ la 6 0) ^sin$cos$ 16 6 1) +sin$ 16 a 0) 

g 
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(4.13) 



where K = fig g/ vf^oS^+^^T^- ^t very early times, just before the interaetion with the pulses, 
the condition 51i/i72 ^ is satisfied, and the dark state \D) -^ \abO), as $(— oo) -^ and K ^^ 1. 
Then, at very late times, the pulses satisfy the condition r22/ili -^ 0, and the dark state \D) -^ \baO), 
as $(+oo) -^ 7r/2 and K ^ 1. 

It is important to note that the strength g of the quantized eavity mode (which is assumed to 
be constant for simplicity) plays an important role here. If the condition g ^ l^o is satisfied, then 
iiT ^ 1 for all times, and the intermediate level \b b 1) is never populated {ilo/g^ 0). This reduces 
the probability of finding a photon roaming inside the cavity, decreasing the chances for cavity decay. 
In addition, if the conditions for adiabatic evolution are fulfilled [TRK92, MYM94], that is 

gT, noT>l and gf7o>r,K, (4.14) 

with T the laser pulse duration, F the radiative decay, and n the cavity decay; then the state vector 
of the system |^(t)) may evolve adiabatically, following closely the dark state \D). 

4.1.4 Density matrix and equations of motion 

For a realistic description of the two-atom + cavity system, in which dissipative channels are ac- 
counted for, we must employ a master equation description. The time evolution of the system is 
described by the following Liouville equation: 

2 

p = -i [H,ff p-p hI^ + ^ Jr, P + Jk/o, (4.15) 

i=i 

where, in the interaction picture and on resonance 

ffeff - -«a^a -'\il l^)^-J-<^l + E (^ l^)j-J-(«l + f \'')n (^l« + H-C-) (4.16) 

is a non-Hermitian Hamiltonian including decay terms from spontaneous emission and cavity decay. 
The supcropcrator Jy describes the return of the electron to the atomic ground states after a 
spontaneous emission, and J^ the corresponding term for the cavity decay [PGCZ95]. 

The decay F terms can be obtained by doing a simple analysis of the spontaneous emission rates 
of the excited levels. For example, the initially populated state |a 6 0) = |a)i|6)2|0)c can only have 
contributions from the radiative decay of states \c b 0) and \a c 0). However, the state \c b 0) is the 
only level which really contributes to the decay, since the state |acO) never happens in this scheme. 
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In addition, this state can also spontaneously emit to other two different states: the state \bbl) and 
some other arbitrary state which is not taken into account in our system evolution. Therefore, the 
contribution made by state \cbO) is r/3. A similar analysis applies to the other possibilities. 

By constructing a density operator with the elements of basis (4.3), plugging it into Eq. (4.15), 
and left- and right-multiplying the resulting equation by the elements of the basis, we find 



Diagonal matrix elements: 



Pabah — ^ [Pabcb — Pcbab) + "IT Pcbcb, (4.17a) 

Pchcb = — ^ — [Pcbab — Pabcb) + -^{Pcbbb — Pbbcb) — ^ Pcbcb i (4.17b) 

Pbhbb — -Z (Pbbcb + Pbbbc — Pcbbb — Pbcbb) + TT \Pchcb + Phcbc) , (4.17c) 

Pbcbc — — ^ — {Pbcba ^ Pbabc) + -^{Pbcbb — Pbbbc) — ^ Pbcbc, (4.17d) 

Pbaba ~ ^ (Pbabc — Pbcba) + TT Pbcbc, (4.17e) 



Off-diagonal matrix elements: 



Pabcb - 2 yP'^'"''' ^ Pabab) - ^ Pabbb - -^ Pabcb, (4.1 (t) 

fll{t) ^ g,^ , ~ \ fA T7 \ 

Pabbb = ^ — Pcbbb + -^ (Pabcb + Pabbc) , (4.1 {g) 

Pabbc ~ — 2 — P"'''"' 2 — ''"''^" ^ 2 ^'''''''' ^ "2 P"-'''"'' *- '' 

f^l(t) ~ ^2it) ^ 

Pabba ~ ^Z Pcbba H ^Z Pabbc, (4-171) 

~ ^i(i)/ . g r ^ 

Pcbab = 2 (P'^'"''' ^ Pcbcb) + 2 Pbbab - ^ Pcbab, (4-17J) 

~ ^iW , Sf \ ^ ~ u ^'7^\ 

Pcbbb — — ^ — Pabbb + -jT (Pbbbb — Pcbcb — Pcbbc) — "jr Pcbbb, (,4.1 IK) 

Qi(t). n^it)^ g,^ . 

Pcbbc = -Z Pabbc H Z Pcbba + X (Pcbbb ^ Pbbbc) — i Pcbbc , (4. i A j 

^i(i) ^2{t) ^g r^ 

Pcbba - — 2 — P"-''''"- 2 — ^"'''"^ 2 ^''''''°' ^ 2 P"''''"-' (4: Aim) 

^iW Sf~ , ~ \ fA -1-7 \ 

Pbbab = — Z — Pbbcb — Z (Pcbab + Pbcab) , (,4.i In) 

Pbbcb — i^ — Pbbab + -Z (Pcbcb + Pbccb — Pbbbb) — -^ Pbbcb, ['i.LIO) 

~ _ ^2{t) gf \ ^~ 

Pbbbc — 7, — Pbbba + Z (Pcbbc + Pbcbc — Pbbbb) — -^ Pbbbc, (4.17p) 

^2(i) ~ g.^ , ^ N I A 1-7 \ 

Pbbba — Z Pbbbc — Z (Pcbba + Pbcba) , (.4.1 Iq) 
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Pbcab = — Pbccb H — Pbaab + ^ Pbbab - -^ Pbcab, (4.17r) 

Pbccb — — -Z — Pbcab -Z — Pbacb + -z[Pbcbb — Pbbcb) — i Pbccb, [iA/S) 

Pbcbb — — ^ — Pbabb + -^ [Pbbbb — Pbccb — Pbcbc) — -^ Pbcbbi l^.i II) 

Pbcba = 2 (Pbaba - Pbcbc) + X Pbbba - ^ Pbcba, (4.17u) 

Pbaab — Z Pbacb Z Pbcab, (4.17v) 

Pbacb 2 — P'""''' 2 — ^'""''' ^ 2 ^''°''''' ^ 2 ^'"^"'^ [-iAiw) 

^2 ft) ~ , g/~ , ~ \ fA -1-7 \ 

Pbabb — 2 P'":!''' + O [Pbacb + Pbabc) , (,4. i IX) 

~ ^2{t) ( . g T ^ 

Pbabc = 2 (/'''cbc - Pbaba) - ^ Pbabb - -^ Pbabc- (4-17y) 

These are the density matrix equations of motion of the two-atom + cavity system, where Ppqrs = 
{pq\p\r s). As before, we have avoided the complex nature of some of these equations by introducing 
the notation 

Ppqrs — iPpqrs, Ppqrs — ^Ppqrs- V^'-Loj 



4.2 Numerical results 

This section presents numerical results obtained by integration of the system of differential equa- 
tions (4.17) using the Runge-Kutta method. Basically, we studied the adiabatic passage method for 
transferring interatomic coherence, when Gaussian and Sech pulses were used to drive a two-atom 
-|- cavity system. For simplicity, we considered a high-Q cavity with constant coupling strength g, 
and no cavity decay (k = 0). For this problem, six different parameters were considered in our 
simulations: the Rabi frequency and the width of the pulses, the time delay, the coupling strength 
of the cavity mode, and the spontaneous emission of the excited states. 

Numerical results corresponding to different values of the Rabi frequency, the pulse width, the 
coupling constant, and the time delay are shown in Table 4.1. 

Table 4.1: Transfer of coherence. Adiabatic passage 



Pulse profile 


F 


^10 


0-1 


r!2o 


(72 


g 


At 


logio P 


Gaussian 


0.00 


2.60 


10.00 


1.00 


1.00 


1.00 


6.50 


-0.68 
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Table 4.1: Transfer of coherence. Adiabatic passage 



14.40 10.10 1.00 4.00 1.00 25.29 -1.54 
0.70 7.00 1.00 2.00 2.00 3.48 -2.19 

0.01 1.10 4.50 1.00 1.00 1.00 0.00 -0.61 
15.00 10.20 1.00 4.00 1.00 25.70 -1.40 
14.80 13.60 1.00 7.00 1.00 36.10 -2.13 

0.10 2.50 7.30 1.00 4.00 2.00 11.69 -1.09 

2.00 6.80 1.00 4.00 4.00 9.80 -1.10 

1.90 6.70 1.00 4.00 6.00 9.40 -1.10 

25.00 30.00 1.00 30.00 1.00 80.50 -1.76 

Sech 0.00 0.90 2.30 1.00 2.00 1.00 0.00 -1.14 

3.80 4.40 1.00 5.00 1.00 11.51 -3.71 

0.01 3.80 4.40 1.00 5.00 1.00 11.48 -2.06 

0.02 5.20 4.10 1.00 5.00 2.00 11.00 -1.85 
6.90 4.10 1.00 5.00 3.00 11.80 -1.85 

As in Ch. 2, we were looking for values of the failure probability of \ogiQp = —2.00 or better. 
During the simulations we fixed the value of 1720 to unity. In this way, we could compare the 
pulses' energies (Rabi frequencies) and widths for the different failure probabilities obtained, and 
then to conclude for which set of parameters we achieved the maximum transfer efficiency. The 
other parameters of the system were allowed to take values from physically reasonable intervals. 
For example, we studied widths in the interval [0, 10] most of the time, going a little bit further in 
some special cases. We also explored time delays in intervals between [0, 50] for Gaussian pulses, 
and [0, 100] for Sech and Lorentzian pulses. 

We first studied the case of no spontaneous emission. We observed from Fig. 4.2 that this 
time there are no sharp dips, but instead we found a long and profound valley consisting of very 
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good transfer of efficiency parameters. It is interesting to note that for Gaussian pulses and no 
spontaneous emission, a good transfer efficiency was achieved for a small value of the pulse amplitude 
compared with its width (see Table 4.1). Fig. 4.3 shows again that Sech pulses are very efficient 
when there is no spontaneous emission. However, bigger pulse amplitudes and longer time delays 
were required, making them less efficient that Gaussians. We tried to extrapolate numerically the 
apparent exponential behavior of these valleys formed by the maximum probability surface, but the 
results we obtained showed no easy connections with the parameters of the system. 

We did not consider Lorentzian pulses in this analysis because they were very computationally 
demanding. 

For spontaneous emission values of F = 0.01 and 0.02, Figs. 4.4, 4.5, and 4.6 showed that Sech 
pulses performed this time better than Gaussian pulses. This was also checked by comparison with 
the ratio of the pulses and the interaction times. This time, the valleys were shallower than those 
for no decay. 

It is very important to note that bigger values of the coupling strength g increased the coherence 
transfer efficiency for no spontaneous emission. However, when the system underwent radiative 
decay, an increment in the magnitude of the cavity mode showed no effect in the failure probability. 
Instead, if the width of the pulses was increased, so was the transfer efficiency. 

In the case of Gaussian pulses, for bigger values of the spontaneous emission, we needed really 
big values of the pulse amplitudes and time delays to get a good transfer efficiency. However, we 
still did not reach the minimum value of the failure probability we were looking for, as shown in 
Table 4.1. 
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^10' ^20 



(a) Surface 



logio{p) 




(b) Contour 

Figure 4.2: Gaussian surface and contour for F = 0.00. Here we have ftio — 0.70, cti = 7.00, 
1720 = 1-00, (72 = 2.00, g = 2.00, At = 3.48, and log^op = -2.19. 
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(a) Surface 



logio{p) 




(b) Contour 

Figure 4.3: Scch surface and contour for F = 0.00. Here we have fiio = 3.80, cti = 4.40, 5I20 = 1-00, 
(72 =5.00, g= 1.00, At= 11.51, andlogioP= -3.71. 
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"lO^"20 



(a) Surface 



logio{p) 




Figure 4.4: Gaussian surface and contour for F — 0.01. Here we have Hiq = 15.10, cti = 13.60, 
^^20 = 1-00, (72 = 7.00, g = 1.00, At = 36.20, and logic p = -2.14. 
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logio{p) 




(b) Contour 

Figure 4.5: Scch surface and contour for F = 0.01. Here we have ftiQ ~ 3.80, ai = 4.40, 5I20 = 1-00, 
a2 =5.00, g= 1.00, Ai= 11.48, andlogioP= -2.06. 
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(b) Contour 

Figure 4.6: Scch surface and contour for F = 0.02. Here we have J7io — 7.30, cti = 5.00, fi20 = 1-00, 
aa =8.00, g=3.00, At ^ 14.70, and logigp = -1.96. 
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4.3 An alternative model 



Consider for simplicity the case in which the two Rabi-frequencies are the same, that is fli{t) — 
il2{t) — ft{t). We assume that the interatomic distance is much less than the wavelength of the 
coherent classical fields, so that the intracavity registration of a photon cannot be used to identify 
which atom is the source of this radiation. The system now has some kind of symmetry and can be 
consider as one consisting of two identical particles interacting with a laser field with Rabi-frequency 
fl(t) and a quantized cavity mode with coupling strength g. The situation is illustrated in Fig. 4.7. 

n(t) 




Figure 4.7: Schematic representation of two qubits interacting with a classical coherent field ft{t) 
and a microcavity mode. The qubits and the cavity can decay with rates T and k, respectively. 

It is interesting to observe that in view of the symmetry of the couplings, the dynamics of the two- 
atom + cavity system can be examined as the evolution of two different systems (one symmetric and 
the other antisymmetric) interacting individually with the classical coherent field and the quantum 
cavity mode, but reciprocally exchanging the interatomic coherence. 

4.3.1 The associated three-level system: EIT 

By adding Eqs. (4.8a) and (4.8e), we obtain 



-t: [Cab + Cba] — ;^ [Ccb + Cbc 



(4.19) 



Then, by adding Eqs. (4.8b) and (4.8d), we have 



-^ [C'cb + Cbc\ — ^-^ [C'cb + Cbc\ H — [Cab + Cba\ + gCbb- 



(4.20) 
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And finally, from Eq. (4.8c), we get 

^a6 = -|[Ceb + Cbc]. (4.21) 

If we assume that the evolution of the system takes place in the Hilbert space spanned by the 
symmetric eigenvectors 

|^i) = -^(|a6 0) + |6aO)), (4.22a) 

|^2) = -^(|c60) + |6c0)), (4.22b) 

\i^3)^\bbl), (4.22c) 

then we define the state vector of the system as the linear superposition 

|*s(t)) = Bi(t)l^i) + B2it)\^J2) + S3(t)|^3). (4.23) 

Here Bi, B2, and B^ represent the probability amplitudes of finding the system in the symmetric 
states IV"!), |'02), and IV's): respectively. These coefficients can be expressed in terms of the individual 
probability amplitudes (4.5) as follows: 

Si(t) = (^i|*s(i)) 

= ^[{abQ\^s{t)) + {baQ\^s{t))] 

= -^ [Cab + Cba] ■ (4.24a) 

Similarly 



B2{t)^^[C,b + Cbc], (4.24b) 

B^{t)=Cbb. (4.24c) 
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On inserting Eqs. (4.24) into Eqs. (4.19), (4.20), and (4.21), wc obtain 



B2 = 
B3 = 






Bo 



r g 

2 2 •" V2 



V2 



So. 



(4.25a) 
(4.25b) 
(4.25c) 



These equations represent the dynamics of a three-level A system with states \ipi) and 1^3) coupled 
to an excited state \tp2) via, respectively, a classical laser field n{t) and a cavity mode field v^g (see 
Fig. 4.8). 

1^2) = ;^(|c6 0) + |&c0)) 




\^,) = j.[\abO) + \baO)) 

Figure 4.8: The associated three-level A system driven by a classical field with Rabi frequency Q(i) 
and interacting with a quantum field of strength y/2g. 

The evolution of system (4.25) is controlled by the non-Hermitian Hamiltonian 



Hs 



which has a dark state given by 



\D) = 



-rj 
n -r V2g 

-^/2g 






(4.26) 



By defining the mixing angle as 



$(t) = tan"^ 



V2g 



(4.27) 



(4.28) 



we can write Eq. (4.27) in the form 



\D) = -sin$(t)|V'i) +cos$(i)|-0; 



3;- 



(4.29) 
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As before, we can transfer population between states Jt/ii) and ji/is) without populating the leaking 
state |i/'2) by employing adiabatic passage via the dark state of the system. The process of coherence 
transfer in this system is illustrated in Fig. 4.9. 

We note that, if \/2g ^ ^, then the state \D) corresponds almost identically to the state IV'i)- 
This means that a single photon excitation is shared among the atoms, favoring the transfer of 
coherence. In addition, the effect of cavity decay is reduced, since the state with a cavity photon 
is not populated. In this limit, a superposition given by the dark state contains only a very small 
component of the single-photon state \b b 1). This increases the lifetime of the combined atom- 
cavity system and is the essential feature of intracavity electromagnetically induced transparency 
(BIT) [FYLOO]. 

4.3.2 The associated two-level system: 27r-pulse coherence transfer 

On the other hand, by subtracting Eq. (4.8e) from (4.8a), we get 

^ [Cab - C\a] = - ^ [Ceb " C\c] I (4.30) 

and by subtracting Eq. (4.8d) from (4.8b), we obtain 

-7: [Ccb — Cbc] ^ —-Z [Ccb — Cbc] — -Z \Cab ^ Cbo\ ■ (4-31) 

Now, if we assume that the evolution of the system takes place in the Hilbert space spanned by the 
antisymmetric eigenvectors 

|^i) = -^[|a6 0)-|6aO)], (4.32a) 

|^2) = ^[|c6 0)-|6c0)], (4.32b) 

then we define the state vector of the system as the linear superposition 

\^A{t)) = Ai(t)l^i) + ^2(t)|^2). (4.33) 

Here A\ and Ai represent the probability amplitudes of finding the system in the antisymmetric 
states \ip\) and |(^2), respectively. As before, these coefficients can be expressed in terms of the 
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Figure 4.9: Transferring of coherence in the associated three-level system. 
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individual probability amplitudes (4.5) as follows: 



A^{t) ^ {^,\^ A{t)) 



V2 



[{ahQ\^A{t))-{haQ\^A{t))] 



— ~~r^ \}~^ab ~ tyfeaj • 



(4.34a) 



Similarly 



A2{t) = -j=[C,t,~Ct,,]. 



On inserting Eqs. (4.34) into Eqs. (4.30) and (4.31), we obtain 



k^-^M. 



(4.34b) 



A. 



2 

0(t) 



Ai A2. 

2 2 



(4.35a) 
(4.35b) 



These equations represent the evolution of a two- level system with ground state \ipi) and excited 
state |(/?2) coupled by a classical laser field il(t), as illustrated in Fig. 4.10. 

I^2) = ^(|c&0)-|6c0)) 



n{t) 



|(^i) = ^(|a6 0)-|6aO)) 
Figure 4.10: The associated two-level atom driven by a classical field with Rabi frequency Vt{t). 

The evolution of system (4.35) is governed by the non-Hermitian Haniiltonian 



Ha = - 



-n 

n -r 



(4.36) 



Here the process of coherence transfer is achieved by a 27r-pulse process, as illustrated in Fig. 4.11. 
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Figure 4.11: Transferring of coherence in the associated two- level system. 
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4.3.3 Transferring the coherence 

By using the symmetric and antisymmetric eigenvectors defined in Eqs. (4.22a) and (4.32a), the 
evolution of the system \abO) — > |6 a 0) can be written in the form 

i=[|^i) + |^i)]^i=[|^i)-|^i)]. (4.37) 

In this way, the transfer of atomic coherence can be viewed as a combination of two different 
processes: an adiabatic passage acting on the symmetric eigenstate \ipi), and a 27r-pulse process 
acting on the antisymmetric eigenstate \(pi). 
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Conclusions 



In this work, wc explored numerically two of the most useful methods for transferring population: tt- 
pulse and adiabatic passage schemes. In particular, we applied these methods to transfer population 
in an atom + cavity system, and to transfer the atomic coherence in a two-atom + cavity system. 
We discussed some important characteristics, as well as advantages and disadvantages of these two 
methods. 

The numerical simulations revealed very interesting resonance-like features in the adiabatic pas- 
sage scheme. By choosing appropriately some parameters of the system, like the ratio of the pulses 
and the pulse delay, we achieved very high transfer efficiencies. To find out different ways of choosing 
the right values for these parameters and reduce the failure probability was the main goal of this 
dissertation. 

It is well known that when a system undergoes no spontaneous emission, we can obtain transfer 
efficiencies of 100% by using 7r-pulses methods. However, when the system undergoes spontaneous 
emission, the transfer efficiency of the method decreases as the radiative decay rate increases. Our 
numerical results showed that efficiencies near to 100% were possible to obtain by using adiabatic 
passage methods for a particular set of parameters, with or without radiative decay. 

For the atom + cavity system and no spontaneous emission, we obtained efficiencies of 99.9% 
transferring the population via adiabatic passage. The Sech pulses showed to be more effective than 
the Gaussian pulses, since they required less relative energy (il/g) and less interaction time {/St/ ag). 
When we consider the effects of spontaneous emission, the transfer efficiency of the 7r-pulses method 
decreased accordingly to the increment in the rate of radiative decay. In particular, for T = 0.10, the 
Gaussian pulses had an efficiency of 93.5%, while the Sech pulses presented a maximum efficiency 
of 92.6%. However, in the adiabatic passage scheme was possible to find a set of parameters for 
which the transfer efficiency achieved a 99.0%. This set was {Vt ~ 2.75, a = 3.09, g = 1.00, Cg = 
2.48, At = 5.29}. Gaussian pulses proved to be more efficient this time than Sech pulses. 

By comparing the ratio of the Gaussian pulses in adiabatic passage, with the ratio of the Gaussian 
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pulses in the 7r-pulses method, and their respective interaction times, we observed that adiabatic 
passage was more efficient transferring population than 7r-pulscs. For Scch pulses to achieve the 
same transfer efficiency, it is required to use more energy and more time delay. 

We also observed that, in adiabatic passage methods, the efficiency can be optimized as much 
as we want by increasing the values of the Rabi frequencies, or equivalently, the widths of the 
pulses. Thus, for bigger Rabi frequencies, bigger efficiencies. Because this is something that has no 
experimental worth, we examined here those values of the Rabi frequency and time delay that can 
be reproduced in a laboratory, and for which very good transfer efficiency values could be obtained. 
So we considered a reasonable value of log^gp = —2.00 as a goal for our numerical simulations. 

We obtained also interesting results transferring the coherence in the two-atom + cavity system. 
With no spontaneous emission, efficiencies of 99.9% were achieved by using the adiabatic passage 
scheme with Gaussian and Sech pulses. However this time, Gaussian pulses required less energy 
and less interaction time than Sech pulses to obtain such transfer efficiency. When spontaneous 
emission was included, Sech showed a better performance than Gaussian pulses with a comparatively 
small ratio of the pulses. A transfer efficiency of 99.1% was achieved for the set of parameters 
{r^io = 3.80, CTi = 4.40, r^zo = 1-00, (72 == 5.00, g= 1.00, At = 11.48}. 

A simple but very useful analytical model used to better understand the transfer efficiencies in 
the adiabatic passage scheme was introduced in Ch. 3. When the widths of the pulses are unequal, 
we do not expect adiabatic passage to work. However, we still found some cases for which a high 
transfer efficiency were obtained. This model described "qualitatively" the dependence of the failure 
probability on the product of parameters fl T. We confirmed that for large values of fi T the 
probability failure power-law decreased. 

Finally, we examined the two-atom + cavity system by using an alternative and simplified model 
based on the superposition of symmetric and antisymmetric eigenstates. We found possible to 
qualitatively explain the transfer of the Zeeman coherence between two atoms, by using an adiabatic 
passage method for the symmetric state, and a 27r-pulses process for the antisymmetric state. 
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Appendix A: Change of variables 



Let us consider the system of three coupled hnear differential equations 

Ca = ^r!i(t)e^"e'^i*C„ (38a) 

d-^r!2(t)e^^e^(^i+^^)*a, (38b) 

C, = 1 f7t(i)e"*"e-*^i* Ca + ]: ni{t)e-'^ e-'''^'+^^^' Cb- (38c) 

It is possible to eliminate the exponentials from Eqs. (38) by making the following change of variables 

Ca^e'^-'Ca. Cb^e'^^^'Cb, C,^e"i'-'C,. (39) 

By differentiating with respect to time, we obtain 

Ca^e'^-'Ca + iKe'^-'Ca, (40a) 

Cb = e'*'* Cb + i^be"^^' Cb, (40b) 

Cc = e">''' Cc + i(j)ce"l'''' Cc. (40c) 

Now plugging Eqs. (39) and (40) into Eqs. (38) gives 

ta = -^c|yaCa + I ^ e-e^(^^+*="^")* Ce, (41a) 






(41b) 



+ . ^IW g^.^^_,(Ai+A.+0.-0,)* ^^_ (4^^) 
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Because we want the phases to cancel out, they must satisfy the conditions 

Ai +(f>,-(f>a = 0, (42a) 

Ai + A2 + 4>c + (jJb -- 0. (42b) 

Solving this system of equations for (j)a, 4>b, and (j)c gives 

(j)a = Ai + (f>c, (pb ^ Ai + A2 + (j)c, 0c arbitrary. (43) 

By choosing 0c — 0, and plugging these phases into Eqs. (41), we get 

iCa{t) = A^Cait) - ^e^" Cc(i), (44a) 

iC'bit) = (Ai + A2)Cb{t) - ^e'P C,{t), (44b) 

^C,{t) = -^e-" CM - ^e-^ C,(t). (44c) 



Finally, the phases a and (i allow us to flip the signs of the terms by choosing suitable values. If 
a = /3 = TT, Eqs. (44) become 

iCa{t) = A^Cait) + ^ Cc(t), (45a) 

iChit) = (Ai + A2)Cb(t) + ^ Cc(i), (45b) 

^Cc(t) = ^ Ca{t) + ^ a{t). (45c) 
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Appendix B: The instantaneous 
Hamiltonian eigenstates 



Let us find the eigenvalues and eigenvectors of tlie instantaneous RWA Hamiltonian (h ~ 1) 



H{t) 









^lit) 








n2{t) 


n^it) 


fi2(t) 






(46) 



Solving the characteristic equation 



det(i7 - Lul) = 



ifli 



1( 



i02 



1^1 i^ 



(47) 



w^- -{ni + ni.)Lu^o 



(48) 



we find the eigenvalues 



--^nl + nl u" = o, 



-^Jnj + nl 



(49) 
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For cj"*" the eigenvector can be determined from the system 



rii z = Qx, 

Q,2Z ^ Qy, 

ilix + il2y — ^ z, 



(50) 



where fl = y^flf + flj. Thus 



X — -prz, y — 7^^7 ^ arbitrary. 



(51) 



Assuming the representation 



, \b) = 



(52) 



the eigenvector associated with w+ can be written as 



!"'"> ^ T2 



>+>>+l=> 



(53a) 



Similarly, for o;*^ and cu , the associated eigenvectors are 



r")^|K-|l^), 



|M^-) 



V2 



' \a) + ^\b} - \c) 



n 



n 



(53b) 
(53c) 



respectively. By using the trigonometric relations 



sin$(t) = -iif , cos$(t) = --^, tan$(t) '^ ^ 



n(i) 



f7(t) 



^2(t)^ 



(54) 
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we can rewrite Eqs. (53) as follows 



\W+) = ^[sin$(t)|a) +cos$(i)|6) + |c)], (55a) 

\W^) = cos$(t)|a) - sin$(i)|6), (55b) 

\W-) = —=[sin^{t)\a) +cos$(i)|6) - |c)]. (55c) 
v2 



81 



Appendix C: The forced harmonic 
oscillator 



The equation of motion for a particle of mass m moving under the combined influence of a linear 
restoring force —kx, a resisting force — 6x, and an external driving force F{t) is given by 

mx + bx + kx = F{t). (56) 

The most general solution to this differential equation is composed of the complementary and par- 
ticular solutions (see [MT95]): 

x{t)^Xc{t)+Xp{t). (57) 



Complementary solution 

The complementary solution has the general form 

Xc(t)-e-^* 



Ci exp i\ Lol - 7^ t +6*2 cxp -i\ cj^ - 7^ i 



(58) 



where 



7=-—, and LUo = \—. (59) 

2m V m 



There are three general cases of interest: 



Underdamping: 7^ < Wg 

Critical damping: 7 — uJq 

Overdamping: 7^ > Wg 

82 



For the case of underdampcd motion, the exponents in the brackets of Eq. (58) are imaginary, and 
the solution can be written as 

Xc{t) = [A sincjit + B coswii] , (60) 



where 



0.1 = ^^0-7^- (61) 

For the case of critical damping, the roots of the auxiliary equation (r^ + 27 r + ujqT = 0) are equal, 
and the complementary function must be written in the form 

Xc{t)^e-''*{A + Bt). (62) 

Finally, for the case of overdamped motion, the exponents in the brackets of Eq. (58) become real 
quantities: 

a;,(t) = e-^*[Cie"i* + C2e-"i*]. (63) 

Particular solution 

Now we seek a particular solution to the inhomogeneous equation in the form 

/•OO 

xp{t)= / F{t')G{t,t')dt', (64) 



where F{t') is the applied external force (inhomogeneity) , and G{t,t') is the Green's function for 
Eq. (56). We define 



^_e-7(t-t')sin[cji(t-t')l ift>t', 

G(i,t')-r""' " (65) 

ifi<i'. 



Then, the particular solution can be expressed as 



xJt)^ f :ffie-T (*-*') sin[cji(t-t')]rft'- (66) 

J-00 "iwi 
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The simple analytical model 

For the simple analytical model studied in Sec. 3.1, we obtained a second-order linear differential 
equation of the form 

^S0[r]+5e[r]=-^'^[T]. (67) 

On comparing Eqs. (56) and (67), we have 

m=l, 6 = 0, fc = l, F{t) = — -$. (68) 

UT 

Upon substitution of these parameters into Eq. (59) and Eq. (61), we get 

cjQ = 1, 7 = 0, LUi — 1. (69) 



For this case we have uJq > 7^ (underdamping). Therefore the most general solution to Eq. (67) can 
be expressed as 



r d$(r') 

[t] = A sin[T] + B cos[t] - / sin[ t-t^] dr. 

J-00 dr 



(70) 
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